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Abstract 

A Trotter-Suzuki mapping is used to calculate the finite-temperature properties 
of the one-dimensional supersymmetric t — J model. This approach allows for the 
exact calculation of various thermodynamical properties by means of the quantum 
transfer matrix (QTM). The free energy and other interesting quantities are ob- 
tained such as the specific heat and compressibility. For the largest eigenvalue of 
the QTM leading to the free energy a set of just two non-linear integral equations 
is presented. These equations are studied analytically and numerically for different 
particle densities and temperatures. The structure of the specific heat is discussed 
in terms of the elementary charge as well as spin excitations. Special emphasis is 
placed on the study of the low-temperature behavior confirming scaling predictions 
by conformal field theory and Luttinger liquid theory. To our knowledge this is the 
first complete investigation of a strongly correlated electron system on a lattice at 
finite temperature. 
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1 Introduction 



Strongly correlated electron systems have attracted considerable interest in recent years 
in view of mechanisms for high-T c superconductivity [I], [|. The t — J model arises in 
various contexts and represents one of the most fundamental systems. The model describes 
the nearest-neighbor hopping of electrons with spin-exchange interaction. The effect of a 
strong repulsive on-site Coulomb interaction is modeled by the restriction of the Hilbert 
space to states without doubly occupied lattice sites. The one-dimensional Hamiltonian 
reads 



+ C j + l,CT C j 



7i,-n i+ i/4), 



where the projector V = Y[Al 



ensures that double occupancies of sites are for- 



bidden. At the supersymmetric point It = J the system was shown to be integrable [|3], [|] 
by the well-known Bethe ansatz |6|]. The ground state and excitation spectrum were 
investigated J7[] and critical exponents calculated by finite-size scaling and conformal field 
theory studies |S], |9| . Here we will discuss the model at finite temperatures. 

According to the seminal work [JTO, 11], |12| thermodynamical properties of general inte- 
grable systems are described by an infinite set of coupled integral equations reflecting the 
existence of infinitely many different rapidity patterns. For the t — J model such integral 
equations were formulated in 



The method consists of a direct evaluation of 
the partition function by taking into account all excited states of the Hamiltonian. The 
excitations are derived on the basis of the so-called string conjecture which describes the 
solutions to the Bethe ansatz equations (for details see In the evaluation of finite- 

size quantities at zero temperature the validity of the string conjecture is controversial 
17] , pll [l9| , ^0| , |2~Tl |22l . Sometimes, a simplistic application is known to lead to erro- 



neous results ||23|| . On the other hand, the string approach to thermodynamics apparently 
yields exact equations. (See for instance J24| where the traditional thermodynamics of the 
quantum RSOS chains were recovered by an independent approach utilizing the fusion 
hierarchy analysis.) For a nice illustration of the successes of the traditional approach to 
thermodynamic properties the reader is referred to the review on the Kondo problem. 

The general problem to be faced within the traditional string approach to thermo- 
dynamics is often the necessity to deal with infinitely many integral equations. In view 
of practical computations this requires a truncation scheme which is difficult to control. 
Also, and more fundamentally, the calculation of quantities other than the free energy, 
for instance correlation lengths at finite temperature, is not possible in the traditional 
approach. 

Our approach takes a different route to overcome these problems. Quite generally, 
<i-dimensional quantum systems at finite temperatures can be mapped onto classical sys- 
tems on d + 1-dimensional inhomogeneous lattices by a Trotter- Suzuki mapping [25]. We 
therefore employ a convenient mapping to a two-dimensional classical model. The clas- 
sical system corresponding to the integrable t — J chain [27] is the Perk-Schultz model 
p8|, |29L BO, which is a well-known multi-component generalization of the six-vertex 
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model. The interesting thermodynamical quantities are expressed by eigenvalues of an ap- 
propriately defined operator, the so-called quantum transfer matrix (QTM). The largest 
eigenvalue, as usual, directly yields the free energy. The further knowledge of the next- 
largest eigenvalues provides the correlation lengths. This represents the important ad- 
vantage [E6l |32|, B3|, |3l, p5f compared to the traditional thermodynamical Bethe ansatz 



requiring all eigenvalues of the Hamiltonian. 

In addition to our general strategy we will incorporate another important ingredient. 
The most common structure of solvable models is the existence of a family of commuting 
matrices comprising the physical operator (Hamiltonian or transfer matrix). We therefore 
utilize a mapping to a classical system where the QTM is apparently embedded in such 
a family. While the general Trotter arguments respect the first conserved quantity only, 
our choice implies the existence of infinitely many conserved quantities. This assures the 
integrability of the QTM from the very beginning. The new approach has been applied to 
several quantum systems |S5|, [H], |23j, [|(], [II] which essentially reduce to scalar 



field theories. In particular a finite set of integral equations has been derived for the spin 
1/2 Heisenberg chain and related models in |3S, fnj. The characteristic feature 

of these equations is the compact formulation originating from the very closely related 
calculation of finite-size eigenvalues of transfer matrices [42]. The investigations of highly 



correlated electron systems, such as the Hubbard model, have just started in j|3| along 
the novel strategy. 

Here we adopt the most sophisticated method [|3^]. The QTM is of the "staggered" 
type and is labeled by two spectral parameters. The dependence on the first parameter 
leads to the commuting family and thereby playing an essential role in the diagonalization 
procedure. The second parameter, on the other side, intertwines the finite temperature 
system with the finite-size geometry of the underlying lattice. This spectral parameter 
represents a spatial anisotropy of Boltzmann weights in the two-dimensional model. Thus 
our procedure can be regarded as a finite temperature extension of the usual Hamilto- 
nian limit. However, we must be careful in taking the limit to obtain finite temperature 
properties: a fine tuning between the parameter and the system size is necessary. This is 
the price to pay for dealing with non-vanishing temperatures. We will overcome the last 
problem by adopting the method developed in the finite-size correction problem ||42|| . All 
information of the Bethe ansatz equations will be transformed into finitely many coupled 
nonlinear integral equations, valid for any system sizes. The Trotter limit then can be 
taken analytically. The resultant equations yield the exact thermodynamical properties 
of the model. 

The study of the thermodynamics of the t — J model using the QTM method has been 
introduced in a short letter j44|. The present paper is devoted to a more comprehensive 
analytical investigation of this system. This work constitutes the first approach allowing 
for analytical and high precision numerical treatments of a multicomponent Luttinger 
liquid. In fact, we confirm the general Luttinger liquid picture of interacting fermions 
analytically in the low-temperature limit. Furthermore, in this limit our equations reduce 
to linear integral equations thus making direct contact with the dressed energy formalism 



for the groundstate investigation of the Hamiltonian El HBL (17 
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The paper is organized as follows. In Section we derive the quantum transfer matrix 
based on the Perk-Schulz model by keeping the integrability structure (see also []48|). In 
Section [$] the eigenvalue equations of the QTM are derived by an algebraic Bethe ansatz. 
In Section |] the eigenvalue equations are transformed into non-linear integral equations 
and general results for physical quantities are presented. Section || deals with several 
limiting cases which are studied analytically. Section |6| contains the conclusion of this 
work. Several more technical, however important aspects of the investigation are deferred 
to two appendices. 



2 Quantum transfer matrix 

The (classical) Perk-Schultz model is considered on a square lattice with periodic bound- 
ary conditions. Each bond of the lattice is occupied by a variable taking on values 1, . . . , q. 
The appropriate non-zero Boltzmann weights TZ 1 ^ associated with a vertex configuration 
a, P, //, v on the lower, upper, left, and right bond 

T^aa( v ) = sinh(r} + e a v)/ sinhry, 

= e a e^smh(v)/smhr], (2) 
Kfflv) = exp(sign(a - /i)v), 



satisfy the Yang-Baxter equation p8| (a, fx = l,...,q). The parameters e a take only 
discrete values ±1 corresponding to bosonic or fermionic statistics of the state The 
Yang-Baxter equation implies the commutation of all row-to-row transfer matrices for 
arbitrary spectral parameters u, v: T(u)T(v) = T(v)T{u) with 



N 



r£(v) = J2U<:T( v )^ ( 3 ) 



/i i=l 



Consequently, the Hamiltonian ([!]) with t — 1 is obtained as the logarithmic derivative at 
v = (where T(0) reduces to the right-shift operator Tr) 



d 



H = —\nT(v 

av 



v=0 

i=l 



£ hi, (4) 



and turns out to be integrable. For the isotropic limit rj (with a rescaling v — > rjv ) 
the non-vanishing matrix elements of the local operators hi read 

yielding the t — J model in the case q = 3 and {ei,e2,e3}= {+ H — }. For q = 3 and 
{ei, €2, £3}= {+ + +} one obtains the Uimin-Sutherland model j|9[ |3| and for q = 4 and 

^2, £3, £4}= {+ H } the Essler-Korepin-Schoutens model [f)0|| , respectively. For more 

details about "generalized t — J systems" see [^l|, ^] . 
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For further algebraic manipulations we introduce the Boltzmann weights TZ and TZ of 
two models related to @ by anticlockwise and clockwise 90° rotations 

According to (§|) we can write 

T(v) = T R e vH+ °^\ T(v) = T L e vH+ °^ 2 \ (5) 

where T is defined in analogy to (|3]) and T RL are the right- and left-shift operators, 
respectively. By means of the substitution 

u = -p/N, (6) 

where f3 denotes the inverse temperature and iV is a large integer "Trotter" number we 
find 

T( U )T(/ /2 = e-WW. (7) 



The partition function of the quantum system 

/ _ \ N/2 

Z = Jim Tr \T(u) T(u)J , (8) 

is identical to the partition function of an inhomogeneous Perk-Schulz model with alter- 
nating rows |48| . The technically more convenient column-to-column transfer matrix of 
such a system is often referred to as the quantum transfer matrix (QTM). Obviously, it 
is a member of the following family of matrices 

N/2 

t qtm (v) = £n<££> + «) <sr( v - »). ( 9 ) 

fi i=l 

at v — 0. The remarkable property of this family of matrices is its commutativity, as the 
TZ and TZ operators possess the same intertwiner which can be derived immediately from 
the Yang-Baxter equation solely for TZ [j52f . 

The free energy / per lattice site is obtained from the largest eigenvalue of T Q ™(0) 

f = ~ lim In ( f j = -\ lim lnA max . (10) 

j3 L^oo \L J (3 n- - 



-•oc 



Furthermore, the next-leading eigenvalues yield the correlation lengths £ of the static 
correlation functions 



- = — lim In 

£ N^oo 



A 



A max 
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Although leaving the evaluation of the finite-temperature correlation length of the t — J 
model as a future problem, we stress that there is, in principle, no obstacle in obtaining 
it in contrast to the traditional string approach. 

Concluding this section we would like to give some comments on the treatment of the 
thermodynamics of Hamiltonian (^) for different particle densities and magnetizations. As 
usual, this is achieved most conveniently by introducing appropriate external fields (chem- 
ical potential /j,, magnetic field h). This in turn modifies @ by an additional factor under 
the trace typically of the kind exp(/3/x £\ n« + (3h ^ s$), where n.j and Si are the particle 
number and spin operators at site i. The definition of the associated quantum transfer 
matrix (§) is modified only by a boundary term |38| depending on the variable hn+i- 



3 Algebraic Bethe ansatz 

As previously mentioned the largest eigenvalue and the next-leading eigenvalues of the 
quantum transfer matrix yield the thermodynamical quantities of the quantum chain. In 
f48fl these eigenvalues were obtained by an application of the algebraic Bethe ansatz. The 



monodromy matrix for this case is defined by an alternating product 



u) ...ni 



(v + u)TZ 



IJ-N A' 



which is related to the quantum transfer matrix (§) by taking the trace over the auxiliary 
space 

>j-QTM(„. \ t,. ri..\ \ " r\i 



(v)=Tr aux £(v) = J2 C > 

A 

The monodromy matrix C(v) satisfies the Yang-Baxter equation 

k${v - «0 4>) ^'H = c » C ^ v ) O - «0- 



(13) 



(14) 



Using the "Neel state" = 1 1, 2, 1, 2, 1, 2, ... ) as a reference state we can construct 
states |\&) by applications of creation operators £f (vi) and £3(101) with different rapidities 
t>i and Wj. These states are shown to be eigenstates |48| of the quantum transfer matrix 
with eigenvalues 



A( v ) = X_(v) + X + (v) + X (v), 



(15) 



where 



A-(«)=n 
a+oo = n 
= n 



w 3 - + i e 1 



v — 


Wj 


v -v k 


-ie 2 


v — 


ffc 


V-Wj 


- ie 3 


v — 





(v — iu — i ei)(t> + i w) 



AT/2 



(u + iw + ie 2 )(u - iu) 



-1 N/2 



(16) 



^+i«) 7v/2 n 



v - v k + i 63 
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and for further convenience we have replaced v by iv. The quantities (3 and jij denote the 
inverse temperature and external fields coupling to the three different quantum states. 
The explicit relation to the chemical potential fi and the external magnetic field h for the 
t — J model {ei, e 2 , £3} ={+ H — } reads 

l_i l = fi + h/2, /i 2 = /i - /i/2, A*a = 0. (17) 

The particular set of spectral parameters {wj, v k } - often referred to as roots or rapidities 
- is determined by the condition that the 'unwanted terms' cancel in the eigenvalue 
problem for This provides an eigenvector with eigenvalue A(v). Using the Yang- 
Baxter equation QUI) it turns out |48| that the parameters {wj,Vk} have to satisfy a 



system of coupled equations. These Bethe ansatz equations read 

Vi + iu + ie 2 \ N/2 yj Vi — Wj - i e 3 y-p «< - v k + i e 3 /3 ( M3 _ /Lt2 ) 



n vi - wj 3 TT 
1 1 ■ — in ■ J- X 



Vi + iu J - L .- L Vi - Wj ±J - Vi-v k -ie 2 

Wj-lU-ieA yr Wj - V k + 1 63 -pr Wj - - 1 63 ^fo^) 

Wi-iu ) \1 Wi- v k J- .J- m;, - Wj + i e x 

It is convenient to use an alternative approach to the Bethe ansatz: The defining relations 
for the rapidities - the Bethe ansatz equations - are equivalent to the analyticity of 
A(v) as a function of v, i.e. the absence of poles in (0). In this sense the denominators 
in (|15D require immediately the Bethe ansatz equations (p~8|) . We like to point out that 



any appropriate treatment of the eigenvalue equations rendering A(v ) analytic may be 
considered as an implicit determination of the unknown roots. This idea is used in the 
next section in a very essential way. 

Next, some general properties of the Bethe ansatz roots are discussed. Consider the 

three different cases {ei, e 2 , e 3 } ={+ H — } and {H h} and { — h +} for which the general 

solutions {wj} and {v k } to ( ]IBD are quite different. The eigenvalues A(v), however, remain 
the same, because each set {e^} describes the same physical system - the t— J model. Thus, 
we may confine ourselves to the case {+ H — } which implies analytical simplifications. 
The largest eigenvalue of the quantum transfer matrix is characterized by N/2 roots for 
each set {wj} and {vk}- Due to the structure of the Bethe ansatz equations (|l^) it is 
natural to assume (at least for vanishing magnetic field h) that the two sets {wj} and 
{vk} for the largest eigenvalue are symmetric with respect to complex conjugation 

Wj = Vj, j = 1, . . . , N/2, h = 0. (19) 

We are led to this conjecture by the typical situation in other Bethe ansatz systems, e.g. 
the groundstate of the Heisenberg model where the roots are symmetrically arranged in 
the complex plane. 

Using (|19D the two different sets of Bethe ansatz equations (with h = 0) reduce to 

v 1 + } u+ 1 \ = yr Vj-Vj + i 
Vi + iu J 1 *■ Vi — Vj 
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Due to the denominator on the right-hand side no real root is allowed, i.e. all roots must 
possess non-vanishing imaginary parts. Eventually, we are interested in the Trotter limit 
N — > oo in order to calculate thermodynamical properties. It turns out that for iV — ► oo 
the roots accumulate at the origin with vanishing imaginary (as well as real) parts. This 
delicate point makes it difficult to analyze the limit N — > oo directly on the basis of the 
Bethe ansatz equations either analytically or numerically. One can overcome this problem 
by introducing well adapted integral equations which is the topic of the next section. A 
detailed investigation of properties of the Bethe ansatz roots is presented in appendix |A]. 



4 Non-linear integral equations 

The eigenvalues A(v) of T Q ™(v) are analytic functions of the spectral parameter v. We 
use this analyticity to determine the largest eigenvalue with explicit representation (|T5| ) 



by a finite set of non-linear integral equations [f£q , |44j| . This approach also allows for taking 
the limit N — > oo analytically. We found that the following combinations of \±(v) and 
X (v) (see (0)) 

A_ 



<B:=l + b, (21) 



A+ + A ' 

b:=y4V' »:=l + 5, (22) 

C - A (A. A + 1l + A ) ' e = 1 + t (23) 

define useful auxiliary functions which satisfy a closed system of functional equations (cf. 
Appendix B). The function b (b) is an analytic complex function along a finite strip in 
the upper (lower) half plane, while c is an analytic complex function on the real axis. The 
analyticity is due to a cancellations of singularities among the A's like in the case of A(v). 
In this sense, the identities between these auxiliary functions encode the information on 
the Bethe ansatz equations as pointed out in the previous section. 

After some lengthy calculations which are described in appendix [B] we obtain the 
following relations for the auxiliary functions 

log b(x) = -2tt/3 ^ b (x + i 7 )+p(ji + h/2) 

- ^ fa * log(l + b)| x . +2i7 - ^ b * log(l + c)| x+i7 , 
log b{x) = -2ixf3 ^(x - i 7 ) + P(ijl - h/2) 

- ^ * log(l + b)U_ 2i7 - ^ * log(l + c)U_i 7 , 
log c{x) = -2tt(3 ty f (x) + 2/3// 

- ^ b * log(l + b)\ x+il - ^* log(l + b)|,_ i7 - ^ c * log(l + c)U 

with the driving terms (and kernels) 

11 2 

2n^ b (x) = - -, 2tt%(x) = and 2vr^ c (x) = ^— , 

x[x — l) X[X + l) X 2 + 1 
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where * denotes the usual convolution f * g\ x = J f{x — y) g(y)dy taken at the indicated 
arguments x, x ± i^y and x ± 2i7 with arbitrary but fixed < 7 < 1. The solution of the 
integral equations provides the largest eigenvalue via: 

log A = -logc(O) + 2/3//. (25) 

For vanishing magnetic field (h = 0) the functions b(x), b(x) are related by complex 
conjugation which leads to a reduction to two nonlinear integral equations. For many 
applications we need only this case to study. 

The integral equations fl24l) can be solved by iteration (for the time being h = 0) 

b( fc+1 ) = exp (-2tt/3 ^ fa + pfi - * b * log(l + W) - ^ fa * log(l + c (fc) 
c (fe+i) = exp ^_2vr/? ^ c + 2/3(1 - 2^ b * log(l + p))) - ^ c * log(l + c (/c) ; 

Choosing appropriate initial functions (see next section) the series {b (fc) ,c (fc) } with 
k — 0,1,2, ... converges rapidly. In practice only a few steps are necessary to reach a 
high-precision result. Moreover, using the known finite fast Fourier transform algorithm 
we can compute the convolutions very efficiently 

In order to calculate derivatives of the thermodynamical potential one can avoid nu- 
merical differentiations by utilizing similar integral equations guaranteeing the same nu- 
merical accuracy as for the free energy. The idea is as follows. Consider the function 

d . . . d . ,, . 1 dc c 

k^-logc with _bg(l + c) = — ^ = — Zc* 

we have 

^ = -^(27^)-^* (Ar c lc 0j ( 26 ) 

which is a linear integral equation for Up. The equations for Ibp and Up (and similar ones 
such as for Zb M , . . . ) are also solvable by iteration after the computation of b and c. 

Now we consider various thermodynamical quantities at intermediate temperatures by 
solving the integral equations (^) numerically. Using the iteration described we obtain 
the grand canonical potential / = f(T,(i) by ( p5|) and its derivative by (p6|), e.g. the 
entropy S, particle density n and the specific heat C are calculated by means of 

As we are interested in the thermodynamical quantities with fixed particle density n we 
have to allow for a temperature dependent chemical potential (i = fi(T) which follows 
from 



-1 

dTj n \dTj\d(iJ T > \dTj n \dT] u \dTj\dn, 



d(i\ _ fdn\ fdn\ fdS\ ( dS\ fdn\ fdn 
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Figure 1: Specific heat as function of T (with h = 0) for different particle densities n with 
n < 0.6 and n > 0.6. 



From the second equation the specific heat for fixed particle density n is obtained within 
the grandcanonical ensemble. For numerical results compare Figure [I]. First of all, we 
note a linear temperature dependence at low T. According to conformal field theory the 
coefficient is given by 7r(l/t> s + l/t> c )/3 where v s and v c are the velocities of the elementary 
spin and charge excitations (see also (|3l|) in section ^1]) . Our numerical data are consistent 
with this expression. (For a completely analytical argument deriving the Luttinger liquid 
properties in the low-temperature limit the reader is referred to section |5T3| ). Furthermore, 
we observe two maxima with changing dominance for increasing particle density n. The 
nature of this structure can be understood from the elementary excitations of the system. 
In the groundstate the particles are bound in singlet pairs with binding energies varying 
from zero to some density dependent value. There are two types of excitations. First, there 
are charge excitations due to energy-momentum transfer onto individual pairs. Second, 
there are excitations due to the breaking of pairs. The latter excitation is of spin type 
at lower excitation energies, but changes the character at higher (density dependent) 
energies to charge type as it describes the motion of single particles. Therefore, the first 
and second maximum at lower densities (figure |].a) are caused by charge excitations due 
to pairs and single particles, respectively. At higher densities (figure |l].b) the maximum 
at lower temperatures is dominated by excitations of pairs whereas the second one at 
higher temperatures is caused by spin excitations. For increasing concentration the spin 
contribution becomes dominant as the charge excitations freeze out. This is in accordance 
with the limiting case n = 1 leading to the spin-1/2 Heisenberg chain. We will come back 
to this point in section 572. The missing spin structure in the specific heat at low and 
intermediate densities is found at quite low temperatures shown in figure |2|. 

It is worthwhile to compare these results with the findings for the Hubbard model 
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Figure 2: Specific heat at low temperatures T for intermediate particle densities n and 
h = 0. 



investigated by the traditional thermodynamical Bethe ansatz |53fl . The structure found 
in the specific heat is explained by spin and charge excitations which do not change their 
character in contrast to the t — J model. For certain densities a low-temperature charge 
peak was found which is caused by single particle excitations. A charge peak at higher 
temperatures appears because of excitations due to doubly occupied lattice sites with 
energies of the order U for large Coulomb interaction. 

To conclude our investigation we present numerical results for other thermodynamical 
quantities. Figure |3] presents the entropy and the compressibility k = dn/dfj, for various 
particle densities. Note the divergent low-temperature compressibility for particle densities 
n — » and n — *■ 1. 




5 Analytical solutions of the integral equations 
5.1 Low-density regime 

In the previous sections we have applied the algebraic Bethe ansatz to the quantum 



transfer matrix and derived non-linear integral equations (|24|) and ( p5|) for the largest 
eigenvalue which is directly related to the free energy fllCf) of the quantum system at 
finite temperature T = 1/ 13. This section is devoted to the low-density limit and analytical 
solutions to the integral equations. 

Let us consider the low-temperature limit for small /i, i.e. the case (3 3> 1 with \i <ti 1. 
Because of 27T/5 ^ t>, c 3> 1 for small values x in (0), we see that the functions b and c are 
almost zero. Therefore, the essential contribution in (fJj|) is caused by x ^> 1. As a first 
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approximation we can assume that b and c do not vary much for which the equations read 
(h = 0) 

log b ~ -13 /x 2 + Py. - log(l + b) - log(l + c), 
log c ~ -2/3 /x 2 + 20n - 2 log(l + b) - log(l + c), 

where the function b becomes real for 7 = 1/2. As a solution we immediately have 

and c= \ , ^-n/^+an , ( 27 ) 



1 _|_ e -/3/x 2 +(3fi l _)_ 2Q-P/z?+Pt* 

which is inserted in (|24| , Pq) again to obtain an approximation for the largest eigenvalue 
(p5|). The analytical expression for the grand canonical potential / = —(log A)//? yields 
(after the change of variable x — > y/P/x) 

--|" 3/2 -6^ r2 + * 3/2 - r2 ^)- < 29 > 

in agreement with the numerical results. This implies for the particle density n = —df/d/j,, 
entropy S = —df/dT and specific heat C = TdS/dT: 

2 _ 2 2 

n = —\/u, S = — T, and C — — T. (30) 
7r 3n on 

These values represent thermodynamical quantities of the t — J model (with h — 0) for 
small particle densities n 1 in the low-temperature limit T 1 such that (T/n) 1. 
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T T 

Figure 4: Comparison of exact (numerical) calculations (denoted by the solid line) with 
the analytical approximation (denoted by the dotted line) according to (pop. Entropy S 
(a) and specific heat C (b) as function of T for different particle densities n. 



The low-temperature behavior of the t — J model described by ( p0|) can be compared 



with independent calculations by means of conformal field theory [pi, 55]. Due to predic- 



tions of conformal field theory we expect for the low-temperature asymptotics 

/ = /o -(^ + pW, (31) 
\6v s 6v c J 

where v s>c and c SjC are the velocities and central charges (c SjC = 1) for the elementary spin 
and charge excitations. For small particle densities we have v S)C = nn [|7j, thus fl29|) and 
fl3"I|) are consistent. 

On the other hand, we can compute the high-temperature limit T — > oo, i.e. /3 — > 0. 
The solutions to the integral equations (p4[) are constants. One obtains 



1 _|_ ' l + 2e^ 

leading to 



and / = -/T 1 log(l + 2e^), (32) 



2(1 — 77,) fi 

S = nlog log(l — n) with fi = T log— r, (33) 

n 2(1 — n) 

as expected by counting the degrees of freedom per lattice site. The simple high- 
temperature limit for finite /i immediately leads to S = In 3 and n = 2/3. 

The analytical benchmarks for T — > and T — > oo imply the consistency of our ap- 
proach. The integral equations for intermediate temperatures have to be treated numeri- 
cally. However, we have found an excellent analytical ansatz for < T < oo valid in the 
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case of small particle densities n C 1 (cf. figure |J). One can extend ( |27| ) to finite ft (with 
7 = 1/2) by 



e -/3/(x 2 +l/4)+j3n e -2l3/(x 2 +l)+2l3v 

= i + e -/3/(^+i/4)+^ and c °( x ) = ! + 2e -/3/(^+i)+^' ^ 34 ^ 

yielding (|25| ) 

= "/I S (^TI73 log + Kix » + ^TT log + c ° w >) ■ (35) 

It turns out that this function describes the thermodynamical properties of the t — J model 
very well. In figure [| the entropy and specific heat (derived from (^)) are depicted as a 
function of T and n. They are compared with the exact quantities calculated numerically. 
In particular, for small particle densities n < 1 the deviation from the exact values is 
negligible. Furthermore, the low-temperature limit provides the integral fl28|) with the 
corresponding asymptotics (p^). 



5.2 High-density regime 

Next, we turn to the study of the t — J model in the limit of large chemical potential \x 
while keeping the temperature finite. With a glance to (|4]) we see that for large \i the 
following scaling behavior sets in 

logb = 0(l), logb = 0(l), logc = /?/i + 0(l). 

In this limit the approximation logc ~ log(l + c) holds up to exponentially small correc- 
tions. Hence, log c can be solved from the last equation in (^4j) in terms of the log(l + b) 
and log(l + b) functions. The actual calculation is done most conveniently in Fourier 
space. For details we refer to Appendix B. The resulting relations are 

log b{x) = -27r/5$(x + i7) +(3h/2 

+ k * log(l + b)\ x - k * log(l + b) |, +2i7 , 

logb(x) = +2vr/3$(x-i7) - fth/2 

+ k* log(l + b)\ x — h* log(l + b)| x ._2i 7 , 

log c(x) = — 2tc(3 k(x) + ft[i 

+ $ * log(l + b)| x _ i7 - $ * log(l + b)\ x+ij , 



(36) 



with the driving term and kernel 

,i k x 



$(x) = , k(x) = — / — 

v ; 2 sinhvnr v ; 2tt 1 H 



+ el fc l 



dk. 
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Finally, the corresponding eigenvalue of the quantum transfer matrix is given by 

, t „, n i f°° log(l + b(x + i 7 )) , i f°° log(l + b(x-i7)) , 
logA = 2/31og2 + /5/i+- / Sl T—r—. I" da?- - / S .\ / . da;. 

2 J-oo smn 7r(x + 1 7) 2 y_ tX) smh7r(x — 17) 



According to this result the particle density in the limit of large chemical potential is 
n — 1 which in fact is the largest possible value. By comparison of ( p6|) with []3S] we find 
the nonlinear integral equations of the isotropic antiferromagnetic spin- 1/2 Heisenberg 
chain. This was expected as the t — J model in the limit of "half-filling" reduces to the 
Heisenberg model. 



5.3 Low-temperature asymptotics 

Let us consider the model at low temperatures but at arbitrary filling. We will employ 
an approximation to obtain the free energy in the low-temperature regime up to explicit 
0(1/(3) terms. The approximation yields a direct relation between the truncated nonlinear 
integral equations and the so-called dressed energy formalism MB, HB[ for the groundstate 



properties of quantum chains. This correspondence itself is novel, and serves as a con- 
sistency check of our results from the nonlinear integral approach on one side and the 
dressed energy formalism and conformal field theory on the other side. 

Consider the system in an external magnetic field. Then b(x) and b(x) are no longer 
complex conjugate to each other. The latter function can be neglected as b ~ 0(e~^ h ) 
for positive h. We choose 7 = 1/2 so that b(x) is a real- valued function on the real axis. 
Then the approximated nonlinear integral equations read 



logb(x) = -(3e° b (x) - ^ fa * log(l + OU+i/a, 

log c(x) = -(3e° c (x) - ^ b * log(l + b)\ x+i/2 - ^ c * log(l + b)\ t 



(37) 



where 



x 



27r* b (z + i/2) - (ji + h/2) and e° t (x) = 2tt^ c (x) - 2/x. 



Numerically, we have verified that b(x) and c(x) exhibit a crossover behavior from 
b(x),c(x) C 1 to b(x),c(x) > 1. Due to the driving term in (|37| ) of order 0((3) the 
crossover becomes very pronounced in the low-temperature limit. We define "Fermi sur- 
faces" by b(±Af,) = c(±A c ) = 1. Keeping these remarks in mind, we split the contribution 
from the first integral term in ([37]) into three pieces 

tfb*log(l + c)U+i/ 2 -> 7T / T SyC (x-x')\ogc(x')dx' 

^ J\x'\>k z 
If 1 

+ — / T s>c (x - x') log(l + — — ) dx' (38) 

27T./ M>At C(x') V ' 

-Ac 



2?r J_ Ac 
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where 

T s>c {x) = 2vr^ fa (x + i/2) = l/{x 2 + 1/4) 



(39) 

T c>c (x) = 2n^ c (x) = 2/(x 2 + 1). 

As the slope is sufficiently steep at low temperatures, we can approximate c(x)(|x| < A c ) 
and l/c(x)(|x| > A c ) by exp(— (log c)'(A c )\x ± A c |) in the vicinity of the "Fermi" surfaces. 
Note that (logc)'(A c ) is of order (3. We can thus justify this linearization over the total 
integral range for the last two integrals in (38). Hence the last term in (p8|) reduces to 

1 r°° 

~— (T SjC (x-A c )+T SjC (x + A c ))x / log(l + e-( lo ^'( A ^)dx 
^ Jo 



{T S)C {x - A c ) + T s , c (x + A e )) 



7T 



(40) 



24(logc)'(A c )' 



The second term in (|38|) can be treated in the same way, and turns out to be identical to 
the result in fllcf ). Similarly, we can linearize the integral equation for c(x). The resultant 
equations are now given by linear integral equations over finite integration intervals 

log b(x) = - (3e° b (x) - (T SjC (x - A c ) + T s , c (x + A c )) 



12(logc)'(A c ) 



— / T. Jx — x') log c(x') dx' 



s,c\ 

\x'\>A c 



2vr 

log c{x) = - pe° c {x) - (T StC (x - A„) + T s , c {x + A b )) 



(T CiC (x-A c ) + T C)C (x + A c ))- 



12(logb)'(A b ) 

7T 



12(logc)'(A t ) 

/ T S:C (x - x') log b(x') dx'-^-f T C:C (x - x') log c(x') dx'. 
27T yu/i>A 6 2tt Ju 1>At 



Let us introduce the short-hand notation 

\ T SiC 1+T C , C J \ logc 
with 

(f> b (x) = -Pe° b {x) - (T s>c (x - A c ) + T s , c (x + A e )) 



(41) 



(42) 



7T 



12(logc)'(A c 



<f> c (x) = -Pe° c (x) - {T S)C {x - A B ) + T s , c (x + A b )) -—-^—— (43) 



(T CiC (x-A c ) + T c , c (x + A c )) 



12(logc)'(A c )' 



We are interested in the evaluation of log c(0). Apparently, it consists of two contributions, 
O(0) and 0(1/(3) which are contained in the integral equation (|4l|) . We will separate these 
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terms by introducing two further 'dressed' functions £5(2;), £c{%) satisfying 

I i?l)*(I)=(t:iV m 



where we have adopted the same abbreviation as in (0). After changes in the order of 
integrations, all integrands are given by products of £&, £ c and known functions. As all 
functions appearing here are even, we arrive at the expression: 



log c(0) ~ 2(/i - 1)(3 + A / e°(x)£„(x) dx + ^ [ e° c (x)Ux) dx 

~ RfWrVfA S ( T ^( A c) -7T- T sA x - A c )f b (a;) dx- — T C)C (x - A,)^(x) dx 

6(logc)'(AJ V 2vr > / N>Ac 2u J\ x \ >Ac 
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6(logb)'(A B ; 



• (r s , c {A b ) -±-J^ T s , c (x - A b )£ c (x) dx) . 



(45) 

By definition (|44}) the contents of the last brackets are nothing but £ C (A C ), £t,(A b ) them- 
selves. We thus obtain the free energy at finite filling as 



f=7T f + ^ [ e° c (x)Ux) dx 

An J\x\>A b Z7r J\x\>Ac 

^b(Ab) tt£c(A 



6/?(logb)'(A 6 ) 6/3(logc)'(A c ) + ° (1//32) - 

Now we compare our results with those found within the dressed energy formalism. 
We quote the result by Kawakami and Yang |4~7| . Immediately seen from their equation 
(2.26), the dressed energy functions for spinon (e s ) and holon (e c ) satisfy the same integral 
equations ( ^2|) except that the right-hand side should be replaced by (— e^, — e°). Therefore, 
we have the connecting relations 

\ogb(x)=(3e s (x) + 0(1/(3), \ogc(x)=Pe c (x)+0(l/P), (47) 

within the 0(1/ (5) approximation. Now that we want to evaluate log b(x), log c(x) and (3f 
including all 0(1/(3) corrections we can replace (log b)'(Ab) and (log c)'(A c ) in the denomi- 
nators of fl46"[) by (3e' s (A b ) and (3e' c (A c ), respectively. Moreover, we find that our £ functions 
are proportional to bulk density functions of spinons and holons: £,b,c( x ) = ^ n Ps,c( x )- T° 



this end compare equation ([44]) with (2.20) in j|7|. Hence the integrands of the first two 



terms in equation (4"6) are products of bare energy functions times bulk density func- 



tions. Therefore, they give "bulk" (or zero temperature) contributions, as expected. The 
remaining terms, which represent finite temperature contributions, can be neatly written 
down by adopting the expressions for the sound velocities of the elementary excitations in 
the dressed energy approach. Namely, let the sound velocities be v a = e' a /27ip a \x a (with 
a = s, c) then the 0(1/ (3 2 ) terms in the free energy read 

7T / 1 1 \ ,, £ 



6(3 2 \v s v c 
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We thus observe various consistency relations. The free energy obeys the prediction by 
conformal invariance and Luttinger liquid theory for arbitrary filling. The sound velocities 
coincide with those obtained from the dressed energy calculations. We remark that there 
should be some subtleties in treating the model for vanishing external fields in the same 
approximation scheme, while the original non-linear integral equations are valid for all 
cases. Even for the linear integral equations obtained in this section analytic solutions 
cannot be found for general filling. Therefore, one has to resort to numerical treatments. 
Here however, we restrict ourselves to the already presented numerical results for the 
original nonlinear integral equations which are not limited to the low-temperature regime. 



6 Discussion 

We have derived eigenvalue equations for the quantum transfer matrix of the t — J model. 
This approach permits the exact analytical as well as numerical calculation of thermody- 
namical quantities. Instead of solving an infinite set of integral equations - as is necessary 
in the traditional thermodynamical Bethe ansatz - we have to solve integral equations for 
only two functions (in the case of vanishing magnetic field). We have considered analyti- 
cally certain low- and high-temperature limits verifying the values predicted by conformal 
field theory. Moreover, the case of intermediate temperatures was treated numerically. As 
shown, the specific heat and compressibility display an interesting behavior in dependence 
on particle density and temperature. Also, the reduction to the Heisenberg model in the 
limit of 'half-filling' has been proved in a straightforward way. A direct relation has been 
established to the dressed energy formulation concerning the groundstate properties. This 
has been easily achieved by a linearization scheme in the vicinity of the "Fermi surface" . 

Let us emphasize again the advantage of the novel method. In the traditional approach 
one has to deal with infinitely many coupled integral equations. Except for special cases, 
the set of equations defies further analysis because of its complexity. Indeed, it was al- 



most 20 years after the derivation of equations |56j that the numerical calculation has 
been done for several thermodynamic quantities of the Hubbard model f5"3|| . Still, the 
explicit calculation allowed for only 2 (!) bound charge rapidities and 15-30 spin rapidi- 
ties, while the original equations contain oo x oo rapidities. Although it was claimed that 
such truncation works for very low temperatures, it might be rather inaccurate at finite 
temperatures unless the Coulomb interaction takes very large values. 

Our treatment, on the other hand, deals with only two integral equations which clearly 
is advantageous in practical calculations. In view of this, the presented work might be 
considered as the first approach yielding explicit and concrete results for a lattice system 
with interacting spin- 1/2 fermions at all temperatures. We actually performed numerical 
computations over a wide range of densities and temperatures. The results are of high 
accuracy (10~ 6 ) in the entire range of parameters. 

The application of the present approach is not limited to the study of ID quantum 
chain problems but also to those of deformed conformal field theories. At the moment the 
main tool in this field is the so-called TBA method which originates from the traditional 
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string approach. As remarked by several authors fl39| , 57|, the novel approach is an efficient 
alternative for scalar models (the sine-Gordon model, spin 1/2 Heisenberg chain). The 
necessity of a generalization of this method to the multicomponent case (Luttinger liquid) 
has been stated as an important problem in f58f . This program has been carried out 
in the present paper. It was shown explicitly for the t — J model. The success of our 
approach is certainly not accidental, i.e. the applicability is not model dependent, but 
rather universal as shown by the general embedding of the QTM into a commuting family 
of matrices. We expect to report in the near future further research on other models, e.g. 
on a generalized Hubbard model []59], as well as on additional thermodynamical quantities 
such as correlation lengths. 
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A Properties of Bethe ansatz roots 

We introduce the function a(x) 

* (v) = — 1 [ ~ ~~t e^. (49) 

\ / J= l J 

According to ( p0|) a root Vk yields 

a{v k ) = -1. (50) 

It is useful to understand the principal behavior of the roots determining the largest 
eigenvalue. This can be done analytically for the case of finite N for fi = h = and large 
(3. Within this limit the integral equation (see appendix [FJ, equations flT5|) and (|76|)) can 
be used to obtain 



loga(vft) = (2A; + l)7ri^$ ( b 7V) (U fc ), 

providing 

v k ~ -~ + 1^1/4 + ^2 + iwcot((2A; + T)ir/N), k = 0, ±1, ±2, . . . (51) 

which will be used as initial values for \i and (3 finite (u = —(3/N). Numerically, these 
values are shown to approximate the actual roots rather well even for finite (3. In the limit 
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N=32, b=8, mu=0 



N=32, b=8, mu=0.5 
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Figure 5: Distribution of roots Vk (full circles) and holes v\ (open circles) in the complex 
plane for iV = 32, (3 = 8, h = and for (a) /i = and (b) fi = 0.5. 



iV — > oo the equation (^1|) reads 

v k ~~ + i v/l/4-i/3/((2fc + l)7r). (52) 

Here one can see easily the accumulation of roots for |fc|^>l. 

In figure [5] the distribution of roots in the complex plane is depicted for two different 
chemical potentials /i. The exact values are computed numerically by the Newton method 
using the initial values fl5l]) . Also shown are the so-called holes v\ satisfying 

a(v h k ) = -1, (53) 

i.e. (|50|), but not coinciding with any roots. As seen from figure |5] the imaginary parts of 
the roots become smaller for increasing [i. Simultaneously, the largest holes converge to 
the real axis. 



20 



B Derivation of the integral equations 

In this appendix we describe the detailed derivation of the integral equations introduced 
in section [|. We define the following auxiliary functions 

A Of 1 , A + + A / c ,x 

a := — , 2l := 1 + o = — , (54) 

A + A + 

b: =X^T' g:=H-^ A - + A + + . + A ° , (55) 

A - A + £ -i | c _ (A- + Aq)(a + + Aq) 

C - Ao(A_ + A + + A )' + A (A_ + A + + A ) ' 1 J 

where Aj are given in flip]). Note that the rescaled eigenvalue (corresponding to zero 
groundstate energy for vanishing chemical potential and magnetic field) 

logA -> log A - N log(l - u), 

(with m = -(3/N) of the QTM is given by 

log A = N log(l + u)-N log(l - u) - log(c(0)) + 2/fy. (57) 

For further convenience we write 

$±(v) = (v ±iu) N /\ q_(v) = l[(v - Wi), q+(v) = l[(v - v t ), (58) 

i i 

which provide the following relations 

X- = q - {V ( + ' ] * + (v)*4v-i), (59) 
q+(v - i) Q_ {v) Q +{v + i)j (60) 

Ao= «-(«+i)^-i) M „ )M „ ) . (61) 

Introducing a, b, 21, 05 as in (|54|-|5"6|) with all subscripts + and — interchanged we find 

A_ + A = A £ <B = A+/b, (62) 
A + + A = A £03 = A_/b, (63) 
(A_ + Ao)(A+ + A ) = A_A + £/c. (64) 

Defining the auxiliary function 

D := — -- (q_(v) $ + {v + i) + g_(t; + i) $+(«)), 

= (A» + Ao(^)) ? 9 '.VI . s = ^fr $ +(^ + i) a ( u )> 
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we can show that D is analytic and non-zero in the upper complex plane C + including 
the real axis. The zeros of the denominator q+(v) cancel with the zeros of 21 because these 
values are identical to the roots V{. Therefore, the only zeros of D are caused by the holes 
in C~ (cf. figure ||, appendix 0). It is useful to rewrite D by means of (|63|) 



D = ^ $+(v)£(v)<B(v), (66) 

= sl±o..( B -i)|tM -L ( 6 7) 

q + {v — i) <P_(f ) b(v) 

which yields ( |64|) 

D3 = $_(!)-i)$ + (t) + i)C/c. (68) 
Now we use the Fourier transform of the logarithmic derivative of f(v) 



f&= I ^ log/(«) e" lto , (69) 
J c ± 2tt L J 

where the integration contour £ is a straight line near the real axis taken in the upper 
and lower half plane such that the contour is parameterized by v = x±i7 with v G C ± , 
x G R and 7 > is calculated for expression fl6"6|) and ©^ for expression (|67D . 

The main idea of this procedure consists in the implicit determination of the unknown 
roots, i.e. the zeros of 21, by employing functional equations in general and the analyticity 
of © in C + in particular. Provided a closed set of equations for b, b, and c in terms of 03, 
05, and £ is found the eigenvalue A of the QTM can be calculated. Employing (|65|), (|66| ) 
and (|67|) we have 



dv 



-1 1 



= e- k q_ + £ + 03, ©« = $ + - q + + £ + 03 = 0, { _ Q) 



2>*>o = e "^- + $+-$- -b, ©« = e fc $_ - e fc g + - b = 0. 

where we used the vanishing of the Fourier transform F(k) for k < (/c > 0) for F(t> ) 
analytic in C + (C - ). Similar relations can be obtained for © reading 



i k>Q = e fc $+-e k q_-b = 0, ©J = e*g+ + $_-$_ 
fCo = $- - 9- + £ + 03 = 0, ©1< = e fc g + + £ + 03. 



(71) 



According to ( |S8|) we have 

® £0 + = e- fc 8 + + € - c, ©go + Co = e fc + I - c. (72) 
Consider the case > 0. With respect to ©•") the equations (|70|) , (|7TD and ( [72]) provide 

e" fc g_ + § = e~ fc $+ - c, (73) 
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(74) 



After inserting the second relation ( |7T|) for it follows 

f = e - k ($ + - $_) - e" fc (» + £), 
c = e - fc ($+ - $_) - e - * (» + £)- <§. 

By similar steps we obtain the full system of equations 

'($+-$_)-(§+£) : fc>0, 

e fc ($_ -$ + ) - e fc (<&+£) : fc<0, 

r k ($+-$_) -e" fc (<B + £) : A;>0, 
($_-$+)-(» + €) : fc<0, 

- fc ($+-$_)- e - fc (25 + £)-§" : £;>0, 
e fe ($_-$+) -e fc (03 + £)-<B : fc<0. 

Furthermore, we obtain: ($+ — $_) = iiVsign(/c) sinh(W). Applying the inverse Fourier 
transform we are led to a system of non-linear integral equations (denoted by the convo- 
lution *) 

logb = - log®- tt b * log £ + /?(// + /i/2), 

log b = <S>f ] - ^ * log 23 - ^ F * log £ + - V 2 ), (75) 

logc = $W -^ fa *log03-^ F *log23-^ c *log£ + 2/3/i, 



with 



= iAM arctan- arctan — J , 



u — 1 f 



u u ) 2irv(v — i) 

v + i f \ _ 1 



= — i AM arctan — arctan — I , \I/ 



M ti 7 b 27Tf (f + i) 



fArl / v — i f + i 

<M j = iAH arctan arctan 



\ m u / 2ir(v 2 + 1) 

The additional terms (3{ji ± /i/2) and 2/3/i in (|75|) are integration constants which follow 
from considering the limit v — > oo 



qPUi-H/2) + I ' " e /3(»+h/2) _)_]_' e /3(M+/i/2) _|_ Qp(jJ.-h/2) _|_ 1 ' 



(77) 



We observe that the Trotter-Suzuki number iV enters only as a simple parameter. In the 
limit iV oo we have 

$ (iV-*oo) = _ 2vr ^^ b - a nd .l.^ 00 ) = -2vr/M c . 



23 



It turns out that the integration contour for \l/ c * log(£ can be moved to the real axis. 
Hence, the equation for b, b is taken on the line v = x ± vy and c on the real axis v = x 
leading to the set of equations (|24]) in section |j. 

Lastly, we want to comment on the treatment of the large [i limit. In this case the 
Fourier coefficients 'c and € coincide up to exponentially small corrections. Hence the last 

equation of ([T4|) can be solved in terms of 23 and 23 

b = sign(A;) — L-j ($ + - $_) + — (23 - S), 
1 + e K 1 + e^ 

S = sign(A0 ($ + " $-) " (»-§), (78) 

c = I = sign(A;) — * (S + - $_) - — Lf (23 + e fc §). 
1 + e |ft| 1 + e^ 

Applying the inverse Fourier transform and respecting the v — > oo limit 



c 



3[i 



K ' V ' ^ e /3/i/2 _|_ e -ph/2 ' ! ' ' ' 



we arrive at the integral equations (|36|). 
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